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Abstract 

Stationary force-free magneto-sphere of an axisymmetric pulsar is shown to have a separatrix 
inclination angle of 77.3°. The electromagnetic field has an Br 1 ! 2 singularity inside the separatrix 
near the light cylinder. A numerical simulation of the magnetosphere which crudely reproduces 
these properties is presented. The numerical results are used to estimate the power of an ax- 
isymmetric pulsar: L = (1 ± 0.1)/U 2 Sl 4 /c 3 . A need for a better numerical simulation is pointed 
out. 
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I. INTRODUCTION 



A neutron star with magnetic dipole p, rotating around its magnetic axis with frequency 
Q loses energy at a about the magneto-dipole rate L ~ p 2 Q 4 /c 3 . This is because the star 
creates free charges that form a magnetosphere with non-zero Poynting flux along open field 
lines - the well-known prediction of Goldreich and Julian [lj. 

The shape of the force-free axisymmetric pulsar magnetosphere was first calculated by 
Contopoulos, Kazanas, and Fendt |2J . This important paper demonstrated that a stationary 
solution does exist, and the power of the pulsar must be close to the magneto-dipole value. 

However, in §2 by solving the stationary force-free equations in the vicinity of the critical 
circle (the intersection of the light cylinder and the equatorial plane), we show that (i) the 
separatrix inclination angle is equal to 77.3°, (ii) the electromagnetic field has an R~ 1 ^ 2 
singularity near the critical circle inside the separatrix. Neither of these properties are seen 
in the numerical results of Q]. 

We therefore repeated the simulation of Q] and found the following (§3). At numerical 
resolution similar to that of Q], our code reproduces most of their results. But at higher 
resolution, the separatrix steepens and fattens, and the singularity of the electromagnetic 
field inside the separatrix starts to develop. The numerical simulation crudely reproduces 
the predicted properties of the magnetosphere. We used our numerical results to estimate 
the power of an axisymmetric pulsar: L = (1 ± 0.1)p 2 Q A /c 3 . 



A. The pulsar magnetosphere equation 

In Appendix A we describe force-free electrodynamics - a remarkable version of plasma 
physics without any plasma properties appearing explicitly. There we also discuss the ap- 
plicability of the force-free approximation. To finish the introduction we must give a brief 
derivation of the pulsar magnetosphere equation jjj. 

It is assumed that electromagnetic forces are much stronger than inertia, thus j x B + 
cpE = 0, or 

V x B x B + V 2 0V0 = (1) 

Here E, B are electric and magnetic fields, p, j are charge and current density, is the 
electrostatic potential; the fields are stationary. For axisymmetric fields, we represent the 
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magnetic field by the toroidal component of the vector potential ip/r and by the quantity 
A = 2I/c, where / is the poloidal current: 



B 



iy <y iy 



(2) 



in cylindrical coordinates r, z; subscripts mean partial derivatives. Using (J2J) in (^Q), one gets 
the pulsar magnetosphere equation |3j: 



(1 - r 2 )(ij rr + -ij r + i) zz ) - -if> r + F(ip) = 0. 



(3) 



here F{ip) = A(ip) , A(ip) is an arbitrary function, = ip follows from the boundary 
condition on the surface of the star, and we use dimensionless units 



c = Q = ft = 1. 

The basic equation Q must be solved with the small distance boundary condition 

„2 



r, z — > 0, 



(4) 



(5) 



( r 2 _|_ z 2^3/2 ' 

which corresponds to the dipole field. 

Michel |4| has found an exact solution of the magnetosphere equation for a magnetic 
monopole rather than a dipole star, proving that solutions of (JB1) which are smooth across 
the light cylinder (r = 1) do exist in some cases. 
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FIG. 1: Stationary force-free axisymmetric pulsar magnetosphere. The thick line shows the sepa- 
ratrix tpQ = 1.27. Thin lines correspond to t/j intervals of O.l^o- 
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II. NEAR THE SINGULAR CIRCLE 



The basic equation Q can be solved in the vicinity of the singular circle, \z\ -C 1 and 
\r — 1| <C 1. In this region, (j3J) can be approximated as 

X{lpxx + ^zz) + ^x = -F, (6) 

where x = r — 1. 

We assume (to be confirmed by numerical simulations of §3) that there is a nonzero return 
current flowing along the separatrix (Fig.l ). Let ipo be the value of the potential on the 
separatrix, and Aq = A(ipo — 0). Then the return current is equal to Aq/2. From (JHJ), we 
get the jump condition across the separatrix 

(WOV^+o - (V^V^-o = ~y x a I ( 7 ) 

In the closed line region, for ip > ipQ, we must therefore have ip ^ (—x) 1 ^ 2 . Thus electric 
and magnetic fields diverge as inverse square root in the vicinity of the singular circle. This 
is an admissible singularity, since the total energy of the fields remains finite. 
We can now find the leading order solution in the closed line region. We set 

ij = ij + R 1/2 f(e), (8) 

where x = Rsin9 and z = RcosO. Then © gives inside the separatrix 

1 d ( df\ 3 , 

' sin ^ +7/ = °- ( 9 ) 



sm9 d9 \ dd) 4" 
Solving this ordinary equation numerically for df/d$ = at 9 = — tt/2, we find f{9 = 
—0.222) = 0. Thus the inclination angle of the separatrix is 77.3°. 

Knowing the separatrix inclination angle, we can solve (JHJ) in the open line region too. We 
assume that in the open line region ^ — ip® = R a f(9), and correspondingly F oc (?/>o~ ^P) 1 
Then equation © reads 

sin 9 d9 \ d9 J sin 9 

Here the two free parameters a and C should be adjusted so as to have (i) /(— 0.222) = 0, 
(ii) /(7r/2) = 0, (iii) no singularity at 9 = 0. Numerical solution gives 

^o-^cxi? 2 - 4 , F(V)oc-(^o-^°- 58 - (11) 

All these properties are roughly reproduced by the numerical simulation presented in §3. 



III. THE AXISYMMETRIC PULSAR MAGNETOSPHERE 



Numerical solution of the axisymmetric pulsar equation Q can be obtained in the fol- 
lowing way . One takes an arbitrary F(ip) and solves (j3J) in the inner (r < 1) and in 
the outer (r > 1) regions using appropriate boundary conditions (the boundary condition 
at the light cylinder being ip r = F(ip)/2). This, of course, does not give the true solution, 
because one gets ^(1 — 0,2:) ^ + 0,2;). One then makes a number of adjustments of 
F(^A, aimed at reducing the jump — 0, z) — + 0, z) for all z. An important finding 
of |2j is that this procedure actually gives an everywhere smooth solution. 

In our simulation, we followed the same method. We were adjusting F(ip) until an 
acceptable solution was obtained. A solution was called acceptable if the following integral 
/ °° dz(ip(l — 0, z) — + 0, z)) 2 1 (1 + z 2 ) was reduced to less then 10~ 7 (starting from « 1 
at F = 0). Equation © was solved by a simple relaxation method. Simultaneously with 
the relaxation, the adjustment of F was carried out, in a way similar to that of Q]. After 
an acceptable solution was obtained, the adjustment of F was stopped, while the relaxation 
was carried out for a sufficient number of steps to ensure that our ip(r, z) does solve (j3J) for 
the obtained F(ip). 

The full function Ff u n which must be used in consists of the regular piece F (the one 
shown in Fig. 2 ), and the delta function piece Fs(ip) = —(J F)5(tp — tpo). In the numerical 
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FIG. 2: The function F which makes the solution of the magnetosphere equation (j2J) smooth across 
the light cylinder. 
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TABLE I: Simulation Results 



Resolution 


(5-function width d 


Separatrix value ipo 


Separatrix inclination, ° 


Power 


200 x 100 


0.02 


1.27 


63 


1.01 


200 x 100 


0.03 


1.27 


60 


1.02 


200 x 100 


0.04 


1.30 


56 


1.06 


200 x 100 


0.06 


1.32 


50 


1.10 


200 x 100 


0.08 


1.35 


48 


1.15 


200 x 100 


0.1 


1.38 


48 


1.19 


100 x 50 


0.1 


1.40 


42 


1.27 


100 x 50 


0.06 


1.34 


48 


1.17 



simulation the delta-function piece was smoothed over the ip interval (ipo,ipo(l + d)). The 
simulation was carried out for different values of d. The figures show the case d = 0.03. The 
table lists some of the simulations that were carried out. The independence of numerical 
results on the radius of the star, for small radii of the star, and on the outer boundaries 
location, for distant boundaries, was checked. 

The power of the pulsar, which is proportional to the spin-down rate, is obtained by 
integrating the Poynting flux over an arbitrary sphere. One gets the dimensionless power 

L = / dMW, (12) 
Jo 

where Aty) = (2 jf dip'F(ip')) is obtained from the regular part of F. From the table, 
we estimate the power of an axisymmetric pulsar L — (1 ± 0.1)/i 2 f2 4 /c 3 . 

Our numerical solution reproduces the features predicted in §2 in the following sense: (i) 
the inclination angle increases with decreasing d and the extrapolated value of inclination is 
about 70° (ii) the maximum of the magnetic field in the inner part of the separatrix becomes 
more pronounced with decreasing d, (iii) the function F demonstrates a singularity similar 
to (fTT|). 

However, we have not accurately reproduced either of the 4 numbers given §2. A really 
good numerical solution should show: (i) the 77° separatrix inclination, (ii) the ip — xpo cx 
R 1 / 2 singularity in the inner region, (iii) the F(ip) oc —(ipo — 4>)°' 58 singularity, (iv) the 
ipo — ij) oc R 2A singularity in the outer region. Until such a solution is obtained, one cannot 



6 



be really sure that a stationary force-free pulsar magnetosphere exists, and one cannot be 
really sure that the power of an axisymmetric pulsar is L — (1 ± 0.1)/i 2 f2 4 /c 3 . 
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APPENDIX A: FORCE-FREE ELECTRODYNAMICS (FFE) 

(This Appendix contains a large excerpt from astro-ph/9902288) FFE is applicable if 



electromagnetic fields are strong enough to produce pairs and baryon contamination is pre- 
vented by strong gravitational fields p. Pulsars, Kerr black holes in external magnetic 
fields, relativistic accretion disks, and gamma-ray bursts are the astrophysical objects whose 
luminosity might come originally in a pure electromagnetic form describable by FFE. 
FFE is classical electrodynamics supplemented by the force-free condition: 

d t B = -V x E, (Al) 

«9 4 E = VxB-j, (A2) 

pE + jxB = 0. (A3) 

V • B = is the initial condition. The speed of light isc=l;p = V- E and j are the charge 
and current densities multiplied by An. The electric field is everywhere perpendicular to the 
magnetic field, E ■ B = 0. The electric field component parallel to the magnetic field should 
vanish because charges are freely available in FFE. It is also assumed that the electric field 
is everywhere weaker than the magnetic field, E 2 < B 2 . Then equation (A3) means that 
it is always possible to find a local reference frame where the field is a pure magnetic field, 
and the current is flowing along this field. FFE is Lorentz invariant. 

Equation (A3) can be written in the form of the Ohm's law. The current perpendicular 
to the local magnetic field can be calculated from equation (A3). The parallel current is 
determined from the condition that electric and magnetic fields remain perpendicular during 
the evolution described by the Maxwell equations (Al), (A2). We thus obtain the following 
non-linear Ohm's law 

. (B- V x B-E- V x E)B + (V-E)Ex B , k . 

J = gr- • (A4) 



Equations (Al), (A2), (A4) form an evolutionary system (initial condition E • B = is 
assumed). It therefore makes sense to study stability of equilibrium electromagnetic fields 
in FFE. One can also study linear waves and their nonlinear interactions in the framework 
of FFE Q. 

One can introduce a formulation of FFE similar to magnetohydrodynamics (MHD); then 
we can use the familiar techniques of MHD to test stability of magnetic configurations. To 
this end, define a field v = ExB/B 2 , which is similar to velocity in MHD. Then E = — v x B 
and equation (Al) becomes the "frozen-in" law 



From v = E x B/i? 2 , and from equations (Al)-(A3), one obtains the momentum equation 



where E = — v x B. Equations (A5), (A6) are the usual MHD equations except that the 



We must mention that the applicability of FFE to pulsars can be questioned If charges 
are not freely available, some regions of the pulsar magnetosphere might exist that should 
be described by a vacuum rather than force-free electrodynamics. While we cannot offer a 
real description of the creation of the space charge, a simple energy estimate shows that the 
system might find a way to put the entire magnetosphere into the force-free regime. Indeed, 
the pulsar luminosity is L ~ B 2 R 6 Q A /c 3 , where B is the magnetic field, R is the radius. The 
number density of charged particles is n ~ QB/(ce). The associated energy density is nmc 2 , 
and the associated power in particles is L p ~ nmc 3 R 2 . The ratio L p /L ~ mc 5 / (eQ 3 B R 4 ) 
is a very small number everywhere in the magnetosphere. Thus, energy-wise, charges are 
indeed freely available. With only a tiny fraction of the pulsar luminosity channeled into 
the charge production, the star will be able to put the entire magnetosphere into the force- 
free state. The real mechanisms for populating the magnetosphere are of course of great 
importance, but these probably include complex interactions of the pulsar radiation, high- 
energy electrons and positrons, the surface of the neutron star, the large-scale and turbulent 
electromagnetic fields — should be difficult to decipher. But FFE might well turn out to be 



9iB = Vx(vx B). 



(A5) 



d t {B\) = VxBxB + VxExE + (V- E)E, 



(A6) 
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a good approximation to reality. 
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